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Quantum diffusion of a relativistic particle in a time-dependent random potential 
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Department of Physics, University of Tehran, Tehran 14395-547, Iran 

We present a rigorous study of quantum diffusion of a relativistic particle subjected to a time- 
dependent random potential with 5 correlation in time. We find that in the asymptotic time limit 
the particle wave packet spreads ballistically in contrast with the nonrelativistic case, which in the 
same situation exhibits superballistic diffusion. The relativistic suppression of wave packet diffusion 
is discussed in connection with statistical conservation laws that follow from relativistic dynamics. 
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I. INTRODUCTION 

Transport in disordered quantum systems exhibits a 
variety of different behaviors ranging from ballistic mo¬ 
tion to diffusion to Anderson localization. In between 
also lies the more exotic type known as anomalous dif¬ 
fusion where the mean squared displacement increases 
slower than a quadratic function of time for ballistic mo¬ 
tion, but not linearly with time as in normal diffusion. 
The other extreme case occurs when the disorder evolves 
in time. Temporal fluctuations of random potential de¬ 
stroy the localization and the transport rate increases 
even beyond the ballistic transport. This effect, which is 
known as stochastic acceleration in classical physics, was 
also predicted in an exactly solvable case of a continuum 
Schrodinger equation where it gives rise to cubic growth 
of mean squared displacement with time [1]. Recently 
this hypertransport has been directly demonstrated in a 
paraxial optical experiment [2]. The cubic behavior can 
also be derived from classical Langevin dynamics with¬ 
out the dissipative term [1] which indicates its classical 
nature. This is expected since the fluctuations of the po¬ 
tential tend to destroy the quantum coherence. In fact, 
the motion of classical particles in a dynamic random po¬ 
tential exhibits this cubic behavior only if the potential is 
5 correlated in time. In a potential with finite correlation 
time the mean squared displacement grows with different 
power ^ = 2.4 at intermediate time (in one dimension) 
[3, 4], and it was predicted in Ref. [4] that the true long 
time asymptotic behavior is normal diffusion. 

Unlike the classical problem spatial correlation of the 
potential also has a considerable effect in quantum dy¬ 
namics. Even in a rapidly fluctuating potential {6 corre¬ 
lated in time) in some cases quantum effects persist for 
long time [5]. This persistence of quantum effects stems 
from statistical conservation laws that govern the quan¬ 
tum dynamics at all times. These conservation laws are 
linked to the statistical properties of the potential. For 
example, the average kinetic energy of a quantum parti¬ 
cle is conserved for certain types of spatial correlations of 
the potential [5] , while for a classical particle it increases 
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linearly with time. As a result of this conservation, the 
diffusion of a quantum particle is suppressed and, instead 
of the cubic growth, ballistic behavior is observed. 

There also exists an earlier analytical study [6, 7] of a 
discrete model, described with a tight binding Hamilto¬ 
nian on a lattice, which leads to normal diffusion. This 
result is due to the microscopic length scale of the lat¬ 
tice which induces a large momentum cutoff. There is no 
such cutoff in the continuum model, therefore the parti¬ 
cle velocity can increase indefinitely. On the other hand, 
indefinite acceleration of the particle in the continuum 
model is the result of nonrelativistic treatment of its mo¬ 
tion. So the other feature that might affect the asymp¬ 
totic behavior is the relativistic limit on the velocity. We 
try to answer the question of how the diffusion of the 
wave packet would be influenced by taking the relativis¬ 
tic effects into account. 

Here we adopt the Dirac relativistic wave equation. 
Besides the fundamental interest, our motivation has also 
been the relevance to various electronic systems which ex¬ 
hibit properties that can be well described by the Dirac 
equation. The main example is graphene, in which elec¬ 
tron transport is essentially governed by the Dirac equa¬ 
tion [8]. Many years after its derivation, the Dirac equa¬ 
tion has attracted enormous attention since the discovery 
of graphene in the last 10 years. Two nonintuitive predic¬ 
tions of the Dirac equation, Kleins paradox and Zitter- 
bewegung, which have been inaccessible experimentally 
can now be tested in a condensed matter experimental 
setup [9, 10]. Even the recent progress in trapped-ion 
experiments has made it possible to simulate such rela¬ 
tivistic quantum effects [11]. Narrow-gap semiconductors 
[12], topological insulators and superconductivity [13] are 
other examples in which the Dirac equation emerges in 
effective description of quantum dynamics. 


II. DIRAC EQUATION AND DENSITY 
MATRIX FORMALISM 

In 1 -|-1 dimensions which we have 1 motional degree of 
freedom the Dirac equation can be expressed using two 
Pauli matrices, ax and cr^: 

d'ip dip 

ih—=-ihcax^+rnc^a„il} + V{x,t)il}, (1) 
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where '0 = ^ j , c is the speed of light, and V(x, t) is assumed to have a zero average, {V{x, t)) = 0, a translationally 
invariant correlation function 


{V{x,t)V{x' ,t')) =v‘^S{t-t')g{x-x'), (2) 

S correlated in time, and an arbitrary correlation in space for which we will later use specific forms for g(x) for explicit 
calculations. Here and below, (■ ■ ■) denotes an ensemble average and {'ijj\ ■ ■ ■ l"^) the quantum mechanical expectation 
value. 

We introduce the density matrix p{x',x,t) = {x' ,t)'tjj{x,t), which is itself a 2 x 2 matrix with elements 

Pij{x',x,t) = 'tjj*{x')'4>j{x), with i,j = 1,2. The time evolution p{x',x,t) = {x' ,t)tjj{x,t) + 'tp'^ {x' ,t)'tfj{x,t) is 

obtained using Eq. (1) as follows: 


pii{x,x,t) 


P22ix',x,t) 


Pl2{x',x,t) 


P2l{x',x,t) 


—c 


—c 






dpi2\ 
dx ) 
dp2i\ 
dx ) 

dp22\ 

) 

dp22 \ 

dx ) 


+ ^[V{x,t) -V{x,t)\pii, 

+ ^[V{x',t) -V{x,t)\p 22 , 

^ - ^ - Pl 2 + — [V{x',t) — V{x,t)]pi 2 

2imc^ 

-^P 21 + -[V{x ,t)- V{x,t)]p 21 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


In terms of new variables p = pn + P 22 , cr = pn — P 22 , t = P 12 + P 21 , 7 = P 12 — P 21 , these equations look simpler; 
also the quantity p is the probability density which is what we need for calculation of moments of position operator: 

p{x', x, t) = -c T + ^-[V{x',t) - V (x, t)]p, (7) 

(T{x',x,t) = -c 7'+ - 4 ^( 2 ^) O]®"! (8) 

f f) f) \ ‘2 i7tic‘^ 1 

T[x',x,t) = -c + - 44(x,t)]r, (9) 

i{x,x,t) = -c a + T + ^[V{x',t) -V{x,t)]a. (10) 


Because the quantities of interest to us are linear in density matrix, we only need the average of density matrix over 
random configurations to calculate the average quantities. To obtain the equations of motion of averaged density 
matrix we need to know averages like {V{y,t)p{x', x,t)). For a Gaussian random variable V{x,t) we can use the 
Novikov’s identity [14]: 


{V{y,t)p{x',x,t)) = J dt" J dx"{V{y,t)V{x\t")) 


Sp{x', X, t) 


^6V{x",t") / ■ 

The functional derivative in the integrand can be calculated first by integrating Eq. (7) with respect to time, 


p{x',x,t) - p{x',x,0) = - f c(^ +-^'jT{x',x,t')dt'+ ^ I [V{x',t')-V{x,t')]p{x',x,t')dt', 


and then 


Sp{x', x, t) 
5V{x",t") 


d d \ 5 t{x', X, t') ! i 


i 

h 


f (— — 

Jt„^\dx dx' 


5V{x",t") 


dt' — / \V{x'^t') — V{x,t 


,^^ 5p{x' ,x,t') , 

SV{x",t") 


( 11 ) 


( 12 ) 


+ - / [5{x' - x") - d[x - x")]5{t' - t")p{x', x, t')dt'. 

^ ./n 


(13) 


In the first two terms, the lower limit of integrals is changed to t" because the density matrix depends only on 
potential at earlier times so the functional derivatives in the integrands vanish for t' < t". The third term is equal to 
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0(t — t")[5(x' — x") — S(x — x")]p(x',x,t), where ff(t) = 1 for i > 0, ^ for t = 0, and 0 for f < 0. Taking the limit 
t" —>■ t, that we need for Eq. (11), and using Eq. (2), we have 


{V{x,t)p{x',x,t)) = - x') - g{Q)){p{x',x,t)), 

{V{x\t)p{x',x,t)) = ^v^{g{0) - g{x - x')){p{x',x,t)). 
By using these quantities we get the following averaged density matrix equations: 


I' 

d 


d d 


-{^{x',x,t)) = -c[ — - — ]{a) + 


dx 

dx\ 

' d 

d 

dx 


d 

d ' 


dx' 

' d 

d 


dx dx' 


-2 7,2 

- x')){"t{x',x,t)). 


(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


With the change of variables X = ^{x + x'), Y = i(x — x'), = ^ + 'W ~ ^ ~ equations are further 

simplified: 


R{X, r, t) = -c^TiX, Y, t) - h{Y)R{X, E, t), 

^(x, r, t) = -C^G(X, r, t) - h{Y)s{x, r, t), 

• c) ^ITTIC^ 

T{X, Y, t) = -c^RiX, Y, t) + —f—GiX, Y, t) - h{Y)T{X, E, t), 

• 0,7 TDC^ 

G{X, E, t) = -c-^SiX, E, t) + —f—T{X, E, t) - h{Y)G{X, E, t), 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where hlY) = ^{g{0) — g{2Y)) and capital Latin letters represent corresponding Greek letters, e.g., R{X,Y,t) = 
{p{x'{X, E), x{X, E), t)). By applying the Laplace transform on t and the Eourier transform on X we obtain 


(s + h{Y))R{K, E, s) = icKf{K, E, s) + Rq, 

— — OliTTiC^ ~ — 

{s + h{Y))T{K, y, s) = icKR{K, Y, s) + G(X, Y, s) + Tq, 

a 

~ f) — OlTflC^ ~ — 

(s + h{Y))G{K, E, s) = -c-^S{K, E, s) + —f—T{K, E, s) + Go, 
(s + h{Y))S{K, E, s) = -cYd{K, E, s) + ^o- 


(24) 

(25) 

(26) 
(27) 


Here Rq, Tq, Gq and Sq are the initial values of the variables (at f = 0) which are known if the initial state of particle 
is determined. Eor every variable we use R{X,Y,s) = e~^^R{X,Y,t)dt and R{K,Y,t) = R{X,Y,t)dX. 

By eliminating variables in favor of G we get the following second order differential equation for G: 


d^G 


5E2 s + h{Y) dY c2 




Yrr?c^ jf? 


(s + h(E))2 + 


G{K,Y,s) = f(K,Y,s), 


where 


^ 1 ^'(^) o , 2mc K{s + h{Y)) - 2im {s + h{Y))'^ 

J{K,Y,s) = -rvTTTT'^o + ^—7— , , , — ttt^Ro - - — , , 


c dY c s + h{Y) 
-\{s + h{Y))Go 


h {s + h{Y))^+ c'^K^ 


h {s + h{Y)Y 


(28) 


(29) 


and h'(Y) = dhiY)/dY. There is one requirement in order to fix the solution of this equation and that comes from 
square integrability of the wave function. Therefore we seek the solutions that at least vanish at infinity. 
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A. Kinematics of wave packet 

So far the problem is reduced to solving the differential equation, Eq. (28). The time evolution of the wave packet 
can be studied through the moments of probability distribution of the position of the particle. Here we restrict 
ourselves to hrst and second moments (multifractal behavior is also expected in the higher moments [4, 5, 15]). The 
moments of probability distribution p can be expressed conveniently as the derivatives of its Fourier transform in the 
following way: 


Mn = ((lA|x"|^/>)) 


J x'^{p{x,x,t))dx 


1 a" 


R{K,Y 


0,t) 


K =0 


(30) 


For n = 0 we simply get the normalization of the wave function at any time {{4’\4’)) = R{K = 0, T = 0, t). Using Eq. 
(24) we have R{K = 0,y = 0,s) = ^o{K-o,Y-0) ^ |.]^g inverse Laplace transform {{iplip)) = Ro(K = 0,Y = 

0) = {^po\1po)^ where tpo is the initial state. Using the formal expression of R in terms of G and initial values, first and 
second moments are obtained as follows: 


QiZinc^ — Q _ \ d — 

Ml = ^y^G(K = 0,y = 0,s) + ^To(K = 0,y = 0)+ i?o(i^,b^ = 0) 


hs^ 

4mc^ d ^ 


K =0 


M2 = 


hs^ dK 


G{K, r = 0, s) 


K =0 


+ ^i?o(if = 0,F = 0)- -^^Ro{K,Y = Q) 


K=0 


(31) 

(32) 


The first term of each of the above equations is unknown and will be determined by solving Eq. (28). We will 
also be using the properties of the Laplace transform to determine the asymptotic behavior of the moments, namely, 
small-s behavior which corresponds to the long time behavior. 


III. EXACT SOLUTION FOR FREE PARTICLE 
WITH INITIAL GAUSSIAN WAVE PACKET 


Evolution of a Gaussian wave packet in the free 

Schrodinger equation is a textbook example. An initial 

i 

Gaussian wave packet, e“U (with zero average 

momentum), spreads during the time evolution but al¬ 
ways remains Gaussian. Calculation of the moments of 
probability density function is therefore simple. For in- 
stance the second moment is given hy p 2 = — + JJiyi t , 
which exhibits ballistic spreading of the wave packet. 
Whereas in general the motion of the wave packet in the 
Dirac equation is more complicated. This is due to the 
existence of two branches of the energy-momentum re¬ 
lation in the Dirac equation corresponding to positive- 
and negative-energy plane wave solutions. Even a initial 
Gaussian wave packet, which in general is a superposi¬ 
tion of both kinds of plane waves, becomes non-Gaussian 
at later times, splits in two parts since positive and nega¬ 
tive components move in opposite direction, and wiggles 
back and forth because of the interference of these two 
parts. These features have already been discussed using 
the numerical solution of the Dirac equation [16] and time 
evolution of the position operator in the Heisenberg pic¬ 
ture [17]. However the time dependence of wave packets 
moving according to the Dirac equation usually cannot 
be determined explicitly. Here we give the closed expres¬ 
sion of the Laplace transform of the density matrix of a 
free Dirac particle with an initial Gaussian wave packet. 

For a free particle (vq —> 0) with a normalized initial 


Gaussian wave packet. 


■00 = fpixR = 0 ) 



a 

p 


e 


l«P + I^P = l, 


(33) 

by inserting the corresponding initial density matrix com¬ 
ponents. 


RoiK,Y)=e-i^'^'e-^^, (34) 

So{K,Y) = (lal^ - (35) 

fo{K,Y) = + (36) 

GoiK,Y) = {a*P - , (37) 

in Eq. (28) and with some algebra, we get 

^ - a^d = {bY + d)e-^'^, (38) 

where , b = ^i\ap - 

\Pp)e~i^^^^, and d = jlfp^(cAr -f 2iRe(a*/3)s) - 
P)s/c^. The general solution for the correspond¬ 
ing homogeneous equation of Eq. (38) is Gie“^-|-G 2 e““^, 
which diverges at both Y = oo and Y = —oo (Gi, G 2 are 
constants). Particular solution has also the same asymp¬ 
totic behavior; therefore it is possible to determine Gi, G 2 
by requiring the complete solution to vanish at infinity. 
















IV. PARTICLE IN A TIME-DEPENDENT 
RANDOM POTENTIAL 
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By doing so we obtain the following solution: 

Xa V2 y' 




16^20 
- {abX^ + Ad) e-“'^erfc 


[abX^ — Ad) e“^erfc [ -I- 


A 


f Xa 

V2y\ 

\2V2 

^ )\ 


(39) 


Now we are able to calculate the moments of the proba¬ 
bility density of the position. Using Eqs. (31) and (32) 
we have 


Ml = 


2cRe(a*/3) y/^Xmc^e 


2mc^ 

h 


s^\/b?s^+Am?e^ 
Re(a*/3)-|-Im(a*/3)s 


xerfc 


A 


4m^c^ 


M2 = 


\2V2c 

2(? 2\f^Xm} 


s'^h^b?s'^-\-Am?d^ 


xerfc 


A 


2^20 


4m 


(40) 


+ (41) 

4s 


To see the long time asymptotic behavior of the moments 
we only need the behavior of the above expressions at 
small s. Keeping the most divergent terms at the limit 
s —>■ 0 and then by the inverse Laplace transform we 

obtain /ii k. 2Re(a*/3) ^1 — v^i 76 rfc( 77 )e’’^ j ct and ^2 = 

^1 — y/Trr]evic{r])e'^^^ where rj = 

It is instructive to see how last two results can also be 
derived from representation of the position operator in 
the Heisenberg picture. Using the Heisenberg equation of 
motion one finds the velocity operator ^ = :^[x,Ho] = 
cax, where Hq = -ihcux-^ + mc^az- For a massive par¬ 
ticle the velocity operator does not commute with the 
Hamiltonian, so unlike the nonrelativistic case the veloc¬ 
ity is not a constant of motion. The equation of motion 
for the velocity can be solved using some operator al¬ 
gebra and then by integrating the velocity one gets the 
following result for the position operator: 

x{t) = x{0)+pc^HQ^t 

+ (42) 

where p = —ihd/dx. The third term is oscillatory in 
time and induces the so-called Zitterbewegung. The mo¬ 
ments now can be obtained as the expectation values of 
different powers of Eq. (42) in the initial wave packet. 
Especially in the long time limit the first moment is dom¬ 
inated by the linear term in time and the second mo¬ 
ment is dominated by the quadratic term as we obtained 
above. The prefactors that we obtained in pi and p 2 
are actually the expectation values {'iI)q\pc^Hq^\'iI)q) and 
(■001 {pc^ ^o" ^) ^ I ^0 ), respectively. 


massless particle. In the case of a massless particle the 
density matrix can be determined exactly too. Because 
even in the presence of a random potential, Eqs. (24) and 
(25) become decoupled from Eqs. (26), (27) and we have 
kK, r, s) = The 

function h{Y) vanishes both at U = 0 and vq = 0; so in 
the calculation of moments of position for which we need 
i? at y = 0 , we get the massless free particle results as 
discussed earlier. However the effect of a random poten¬ 
tial can be observed in correlations of the wave function 
at two different points (Y 7 ^ 0). Further with the inverse 
Laplace and Fourier transforms we obtain 

R{X,Y,t) = ^^e-'‘(^)*([I-2Re(a*/3)]e"^(^+^*)' 

V 27rA 

-f[I-b2Re(a*^)]e"^(^"^*)"). (43) 

We can see that at later times the wave packet decom¬ 
poses into two smaller wave packets moving in opposite 
directions at the speed of light. This is actually what 
happens in the free Dirac equation. Here in addition we 
have an exponential factor, due to the random poten¬ 
tial, representing the decay of correlations of the wave 
function at two points at the distance Y with the time 
scale To = 77 (Yj which depends on the spatial correlation 
function. We may also note that this characteristic time 
diverges both at U = 0 and at the zero disorder limit 
Vo = 0 . 

massive particle. For this general case the differential 
equation, Eq. (28), is rather complicated to be solved an¬ 
alytically but the numerical solution is possible, although 
imposing the boundary conditions is not straight forward 
numerically. As we explained above the solution must 
vanish at infinity but numerically we can only determine 
the solution at a finite interval. At the boundaries we 



Y 

FIG. 1: (Color online) Numerical solution of Eq. (28) for 
r = 1 , A = 0 , s = 0 . 001 , a = ~ \[\ other 

parameters equal to 1. The real part is odd therefore vanishes 
at y = 0 and the solution is independent of s at small s. 
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FIG. 2: G'{Y,s) for r = 1, s = 0.01 and other parameters 

equal to 1. Inset is the plot of —G'{0, s) in log-log scale which 
shows the i divergence in s —>■ 0. 


G\Y,s) 


dK 


G{K,Y,s) 


K=0 


and the inset is a log-log 


plot of —G'{Y = 0, s) which indicates power law singu¬ 
larity at s =_0. By fitting the inset to s~^ we obtain 
(5 = 1, i.e., G'{Y = 0,s —0) = —£2^; therefore the 
first term in Eq. (32) has also s~^ singularity at s = 0 
therefore at t —>■ oo 


fJ-2 



^ lim sG'{Y 
h s^o 



(44) 


For comparison the corresponding nonrelativistic result 
(Eq. (17) of Ref. [1]), with parameters that used here for 

r = 1, is /X 2 = ^ -h We repeated the 

calculations for other values of r > 1 and obtained the 
same exponent i5 = 1 and the same asymptotic behavior 
of moments. 


impose (we use MATHEMATICA for numerical solution of 
differential equation) a very small value for G but we need 
to make sure that the interval is large enough so that the 
result does not depend on the boundary values. First we 
note that in order to calculate the moments (31) and (32) 

we only need G at iY = 0 and ^ at K = 0, respectively; 
so first we obtain the corresponding equations of these 
quantities from Eq. (28), which will be relatively simpler. 
Moreover we are interested in the long time behavior of 
moments; therefore we need to know the s dependence of 
the Laplace transforms at s —>■ 0. We also take the short 
range spatial correlation g{x) = , where ^ is the 

correlation length and r is a positive integer. 

Figure 1 shows real and imaginary parts of G{K = 
0,Y,s) for r = 1, s = 0.01, a = and 

other parameters equal to 1. The real part vanishes 
at Y = 0, which is necessary because otherwise we 
obtain a complex value for /ii [see Eq. (31)]. More¬ 
over the solution converges as s goes to 0; so from 
Eq. (31) and the inverse Laplace transform we have 

/j,i = c ^1 — ^^^ImG(iY = 0,Y = 0)^ t. Figure 2 shows 


V. SUMMARY AND DISCUSSION 

To summarize, we have studied the evolution of a wave 
packet in the Dirac equation using the density matrix 
formalism. We derived exactly the density matrix for 
a massless particle in a time-dependent random poten¬ 
tial as well as for a free particle with an initial Gaussian 
wave packet. For a massive particle in a random time de¬ 
pendent potential we obtained the asymptotic behavior 
of average position and mean squared displacement. In 
contrast to superballistic diffusion of a Schrodinger par¬ 
ticle for the correlation function g(x) = with 

r = 1 [see Eq. (2)], our results exhibit ballistic motion 
in this case [Eq. (44)] . We find ballistic behavior in long 
time limit for r > 1 as well, which coincides with nonrel¬ 
ativistic results. 

The difference between two cases suggests the existence 
of another conserved quantity in the relativistic case. A 
nonrelativistic quantum particle gains energy from the 
Gaussian correlated potential (r = 1) as its average ki¬ 
netic energy is not a constant of motion. For a relativistic 
particle we have 


(Ek) = ((moW) 


J dxitp^ {x)Ho'il^{x)), 


dx ^ {-tp*{x){-ihcaxdx + mc^(Jz)ij'ipj{x)), 
ij'=l 



-ihc—{T{x', X, t)) + m(?{a{x', x, t)) 




-^-^nX,Y,t)+mc^S{X,Y,t) 


(45) 

(46) 

(47) 

(48) 

(49) 
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In going from Eq. (47) to Eq. (48) we have performed the change of variables X = + x'),Y = ^{x — x'). In 

the third line, the integral of the term -^T{X,Y,t) vanishes because T{X = ±oo,Y,t) = 0. Also note that in the 
integrands the differentiations are being done first and then the values x' = x {Y = 0) are set. Now by taking the 
time derivative and substitution from Eqs. (22), (23) and using h{fi) = h'{0) = 0, we obtain 




ihc d 


f{X,Y,t)+mc^S(X,Y,t) 


r=o 



'iHc^ d^R 

2 dXdY 


,dG ihc d 

W 


{h{Y)T) 


itu? dR 

~1~W 


X=oc,Y=0 


riC 

—mc^— - mc^h(Y)S 

oY 


Y=0 


ihc^ dR 


2 dY 


X=-oo,Y=0 


= 0 , 


(50) 


(51) 

(52) 

(53) 


which shows that the average kinetic energy is conserved 
even for the Gaussian correlation function. The explicit 
derivation of density matrix that is obtained in the mass¬ 
less case (Eq. 43) reveals an interesting feature. The cor¬ 
relation function {'tjj''^{x',t)ip{x,t)), in addition to tran¬ 
sient time dependence due to the motion of the wave 
front, exhibits an exponential decay with a time scale 
To = / ,N induced by spatial correlations of a time 


dependent random potential. 

VI. AKNOWLEDGEMENT 

I would like to acknowledge the hospitality of ICTP 
where part of this work was completed. 


[1] A. Jayannavar and N. Kumar, Physical Review Letters 
48, 553 (1982). 

[2] L. Levi, Y. Krivolapov, S. Fishman, and M. Segev, Na¬ 
ture Physics 8, 912 (2012). 

[3] L. Golubovic, S. Feng, and F.-A. Zeng, Physical review 
letters 67, 2115 (1991). 

[4] N. Lebedev, P. Maass, and S. Feng, Physical review let¬ 
ters 74, 1895 (1995). 

[5] A. Jayannavar, Physical Review E 48, 837 (1993). 

[6] A. Madhukar and W. Post, Physical Review Letters 39, 
1424 (1977). 

171 A. Madhukar and W. Post, Physical Review Letters 40, 
70 (1978). 

[8] K. Novoselov, A. K. Geim, S. Morozov, D. Jiang, M. K. 1. 
Grigorieva, S. Dubonos, and A. Firsov, nature 438, 197 
(2005). 

[9] M. Katsnelson, K. Novoselov, and A. Geim, Nature 


Physics 2, 620 (2006). 

[10] M. Katsnelson, The European Physical Journal B- 
Gondensed Matter and Gomplex Systems 51, 157 (2006). 

[11] R. Gerritsma, G. Kirchmair, F. Zaehringer, E. Solano, 
R. Blatt, and G. Roos, Nature 463, 68 (2010). 

[12] W. Zawadzki, Physical Review B 72, 085217 (2005). 

[13] O. Vafek and A. Vishwanath, Annual Review of Gon- 
densed Matter Physics 5, 83 (2014). 

[14] E. Novikov, Sov. Phys. JETP 20, 1290 (1965). 

[15] J. Bouchaud, D. Touati, and D. Sornette, Physical review 
letters 68, 1787 (1992). 

[16] B. Thaller, The dirac equation, vol. 31 (Springer-Verlag 
Berlin, 1992). 

[17] A. Barut and A. Bracken, Physical Review D 23, 2454 
(1981). 














